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Abstract 


A computational procedure for describing transport of electrons in 
condensed media has been formulated for application to effects and 
exposures from spectral distributions typical of electrons trapped in 
planetary magnetic fields. The procedure is based on earlier param- 
eterizations established from numerous electron beam experiments. New 
parameterizations have been derived that logically extend the domain of 
application to low molecular weight (high hydrogen content) materials 
and higher energies (~50 MeV). The production and transport of high 
energy photons (bremsstrahlung) generated in the electron transport 
processes have also been modeled using tabulated values of photon pro- 
duction cross sections. A primary purpose for developing the procedure 
has been to provide a means for rapidly performing numerous repetitive 
calculations essential for electron radiation exposure assessments for 
complex space structures. Several favorable comparisons have been 
made with previous calculations for typical space environment spectra, 
which have indicated that accuracy has not been substantially compro- 
mised at the expense of computational speed. 

Introduction 

The constituents of ionizing radiation fields in the space environment of most importance for effects 
on human or other biological systems and radiosensitive electronic components are the high-energy nuclei 
of the galactic cosmic rays (GCR), the energetic light ions of sporadic solar particle events (SPE), and the 
protons and electrons trapped within magnetic fields associated with a stellar or planetary object. Inter- 
actions of these radiations with condensed matter can produce secondary fields of high-energy neutrons 
and photons. Prediction of the nature of effects caused by such radiations depends upon understanding 
the manner in which the radiation fields propagate through various materials comprising the media of 
interest. A comprehensive exposition of various aspects of light and heavy ion transport is given by 
Wilson et al. (ref. 1) that emphasizes the adequacy of a one-dimensional approximation in the attenuation 
of such ions in condensed matter. The error is proportional to the second power of beam divergence to 
radius of curvature of the shield. Mainly, low-energy neutron transport requires added consideration. 
The much lighter electrons, however, are much more subject to undergo change in their direction of 
motion when they interact with the charged particles comprising the atoms and molecules of a substance. 

An excellent treatment of the various complexities of detailed electron transport is given by Jenkins, 
Nelson, and Rindi (ref. 2), which is devoted to the analytical description of electron transport by Monte 
Carlo techniques. Such statistical approaches have the advantage that the physical processes involved 
may be modeled in detail, and comparisons with controlled experiments are usually very favorable so that 
when the physical models have been verified, Monte Carlo results are often used as “benchmark” calcu- 
lations. Some disadvantages of the Monte Carlo approaches are that they often tend to require substantial 
computer storage and execution time. This is most important in an engineering design environment 
where optimization processes depend on high-performance computational methods. 

The present work describes a formulation in which the philosophy of one-dimensional transport is 
retained, but modified to try to account for average effects of multiple scattering events and high fre- 
quency energy fluctuations. Complex mathematical descriptions of individual collision processes have 
been replaced with a continuous, deterministic representation. Thus, the computational complexity and 



mass storage requirements of Monte Carlo techniques are greatly reduced. Nevertheless, as with most 
engineering approximation approaches, some aspects of the physical phenomena are not represented, and 
the range of validity of the approximations invoked should be carefully defined. The purposes of the 
present study are (1) to describe the formulation of the computational procedure in some detail, along 
with input data base structure; (2) to compare some results with suitable benchmark calculations and 
define some limitations and ranges of validity; and (3) to demonstrate by means of sample calculations 
some representative applications of the procedure. 

Electron Transport 

Computational Formulation 

The paths of sufficiently energetic electrons traveling in condensed media may be roughly approxi- 
mated by straight-line motion. In the absence of catastrophic interactions (e.g., positron annihilation, 
nuclear reaction processes), the straight-ahead approximation becomes more favorable with increasing 
energy. The electron transit through a substance may be characterized by the stopping power -( dEldx ), 
or loss of energy E, per unit distance x. Approximate stopping power formulas have been developed from 
first principles, and later extensively refined to incorporate corrections to earlier expressions. 
An excellent overview on stopping power theory and application in the section by Berger appears in 
reference 2. 

If an electron having initial energy E is assumed to lose its energy continuously in transit, its maxi- 
mum travel distance or path length R may be found as 



1 

(=&L\ 

\ dx ] 


dE' 


(1) 


and is referred to as the Continuous-Slowing-Down-Approximation (CSDA) range, or pathlength. In 
practice, straggling, and especially multiple scattering, makes the pathlength a complicated set of 
randomly directed line segments. The actual penetration depth or range of electrons is observed to be less 
than the CSDA range, with the deficit dependent on both material composition and energy. Tabata, Ito, 
and Okabe (ref. 3) have developed a parameterization for the practical range of electrons based on 
experimental tests involving several materials and a broad energy range (0.3 keV to 30 MeV). Tabata 
explicitly parameterizes the practical (extrapolated) range R sx of an electron in material of atomic number 
Z and atomic weight A as 
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where x = El(m e c - ) with m e as electron rest mass and c the speed of light. The coefficients a ,• are given in 
terms of A, Z, and constants bj= i _9 as 
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With the bj specified as 
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the formulas may be used to calculate R sx in units of g/cm 2 . 

The parameterization of reference 3 has been examined as to the relationship with the CSDA range 
computed from the stopping power tabulations of Berger and Seltzer (ref. 4). Plots of the CSDA stopping 
powers and corresponding range values are given in figures 1 and 2, respectively, for selected elemental 
species and for electron energies between 0.01 and 100 MeV. When the parameterized practical ranges 
are ratioed to the CSDA ranges, values generally less than unity occur over the quoted energy range 
of applicability for the middle and higher atomic number elements (Z > =10), as expected. However, 
the formulas of reference 3 lead to ratios exceeding unity for these heavier elements at higher energies 
(> =20 MeV). For lower Z elements, the CSDA range is exceeded at even lower energies. Since the for- 
mulations of reference 3 are based primarily on data citations for the heavier elements, it was felt that a 
reevaluation be considered using the observed range ratios for the heavier elements along with a rational 
extrapolation to higher energies and lower atomic weights. 

It is convenient to define the ratio Fr = R ex /Rc sda- Plots of the semiempirical functions of refer- 
ence 3 used in CSDA range ratios suggest that lower order polynomials are suitable functional forms for 
elements with Z> =13. We have parameterized Fr as cubic equations of the form 

F R = a(Z)X 3 +b(Z)X 2 +c(Z)X+d(Z)m R<iX (3) 
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in terms of a normalized energy variable, 


x= l°£l0 (^Mev)~ i°£l0 (0-0l) 

= [log 10 (100)- log 10 (0.01)] 

It is found that the coefficients exhibit approximately regular variations with Z and may also be repre- 
sented by simple mathematical functions. Reasonable fits to the coefficients (including rational extension 


to low atomic numbers) have been found as: 

a = 0.02258(y z f -0.1 3854(y z ) + 0.3695 (5) 

b = -0.00605(.r z ) 3 + 0.04388(.x z ) 2 + 0.07536(x z ) + 0.075 (6) 

c = — 0.56 Jl .0 - exp -0.0030706(y z ) 41514 J (7) 

d= 0.971 -0.07l(y z ) (8) 


with yz = Zh 2 and x Z = y z - 1 . 
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In order that the new parameterizations Fr are constrained to be less than or very nearly equal to 1 , 
they are modified by an exponential “filter” such that when Fr approaches 1 at higher energies, the poly- 
nomial is suppressed exponentially and replaced with a corresponding approach to unity. The exponential 
function is of the form 


qp(X) = exp -A(xy 


( 9 ) 


where X is the normalized energy variable previously defined. For the present application, the constants 
A and n are chosen so that the slope of the function is -80° at tp(X) = 0.5. These constants may then be 
expressed in terms of a “transition energy variable” Xj as 

« = 16.364X r 
A = 0.693 15(X r )“" 

The transition energies for which these conditions occur are somewhat arbitrary, and have here been 
selected as a function of atomic number so as to produce a reasonable approach to the CSDA range with 
increasing energy. The functional form chosen is 

X T {Z) = -0.0043438X| + 0.078366X Z + 0.5753 (11) 

where Xz = Z^ 2 - 1. Thus, the final form of the expression for practical (extrapolated) range as used in 
the present study is 



^x/^SDA = Fr{X, z) cp(x Z) + [l - cp(X, Z)\ ( 1 2) 

Graphical representations of the range ratios for the earlier and present formulations are given in 
figure 3. In this way, the extrapolated range approaches the electron pathlength at high energies, and, for 
light elemental targets, reaches the asymptote at lower energies. 

In order to account for the differences between observed energy dissipation by electrons and that pre- 
dicted by CSDA, Kobetich and Katz (refs. 5 and 6) have devised a relationship of the form 

G{E j).^n (i3, 

dt 


where G is the energy dissipation of an electron of initial energy E having residual energy W, after 
traveling distance f; r| is the probability of transmission at distance t. Their parameterization of r| in terms 
of range R is given as 
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and is limited to values of Z > 1 . 
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In the present formulation, the formula for p of Kobetich and Katz has been modified as 

p = 1.2 + 1.3415 exp (-0.04723Z) (16) 

to extend the parameterization logically into the neighborhood of Z = 1 (where the original formula tends 
to become indeterminate) and to retain the character of the formula at higher Z values. It was found that 
this modification appeared to be a viable representation. The effects of path-length straggling near the 
distance at which the electron has lost most of its energy are also parameterized in reference 5 to constrain 
the dissipation function G as 


G{E,t>0.9R ) 


G(E,t = 0.9R) 


r]{t) 
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(17) 


This expression has been retained in the present formulation. 

For a fluence (or flux) of electrons differential in energy, fy(E,l=0 ) — having initial energies E, the dose 
D, at depth t — is given in terms of the dissipation function (ref. 7) as 

D(t) = (1/ 2 k) f G(E, t)<p(E, t)dE (18) 

The dose may also be written in terms of the residual energies W as 

D(t) = f S(W)<b(W,t)dW (19) 

(Note that the energy dissipation function has been derived on the basis of electrons incident on a semi- 
infinite slab for total solid angle of 2 tc sr.) Equating the two integrands and dividing each side by 
differential thickness element dt leads to an expression for the electron spectrum at location t in terms of 
the initial spectrum: 


<KM = (1/2jc) c|)(£, 0)G{E,t) f E ) T (20) 

where S represents the CSDA stopping power. The residual energy W at distance t, for an electron of 
initial energy E and having range R, is found by solving the equation 

R ex (fV) = R ex (E)-t (21) 

for W when the practical range satisfies the condition R QX (.W) > 0. This operation is performed 
numerically in the present algorithm. 

As electrons undergo deceleration in transit through condensed matter, energetic photons 
(bremsstrahlung) are produced that also propagate through the medium and contribute to the general en- 
ergy deposition (dose). The photons are generally distributed spectrally from very low energies up to the 
electron energy producing them, and their spectra are related to production cross sections that are depend- 
ent on both the atomic composition of the material and the electron energy. When electrons of energy W 
at distance t produce photons of energy in the interval between E y and E v + dE v , the source term t, (in, 
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for example, units of photons/cm 3 -MeV), at distance t in the material, may be written in terms of the 
differential photon production cross section, o(W, E v ) 

‘Q(E v ,W) = jj (0 (P (W',t) o(W’,E v ) dW (22) 

where W' > E v . The cross sections o(W, f? v ) are generally complicated functions of W and E v , and 
material composition. They have been extensively tabulated by Seltzer and Berger (ref. 8) for wide 
energy ranges and most elements of the periodic chart. The effective production cross sections for a 
given material are determined in the present calculations by appropriate spline interpolations of 
the Berger-Seltzer tabulations. Plots of scaled cross sections for selected atomic species are given in 
figure 4 ((a)-(f)). The differential cross sections have been scaled with multiplication by the factor 
i? v (|3/Z) 2 , with p being the ratio of electron speed to the speed of light. Individual curves represent cross 
section variation for a constant ratio of emitted photon energy to incident electron kinetic energy. 

The attenuation of the emitted photons in the material can be characterized by a total extinction coeffi- 
cient pr and the photon differential energy spectrum (p v at distance t in the medium may be found by 
using the radiative transfer equation: 

<t> v (£v,t) = J XiX’Ev) exp [-p r (/-x)]tft ( t>x ) (23) 

Photon energy deposition contributing to the general ionizing radiation dose may be characterized by an 
“energy deposition” absorption coefficient p£ in the calculation of dose due to bremsstrahlung, 

D v(t) = f VE E v <M £ W) dE v (24) 

The energy absorption coefficient is generally less than that for total attenuation pj. Photon absorption 
coefficient data used in the present procedure are shown for selected species in figures 5 and 6, and have 
been taken from the data tabulations of Storm and Israel (ref. 9). The present code formulation assumes 
all photons generated propagate in the direction of electron motion. This approximation should result in 
overestimates, or conservative values, for photon energy deposition in a given direction. 

Comparisons With Other Calculations 

The present computational procedure was used to replicate several calculations performed with other 
electron transport codes. In what was considered to represent a benchmark comparison, the Monte Carlo 
code TIGER-P (ref. 10) was used to generate dose versus depth curves in water for electrons and associ- 
ated bremsstrahlung. An initial spectrum based on the 10-day mission fluence of STS-63 was propagated 
into the EGO medium at normal incidence. The results of this generally favorable comparison, along with 
the initial spectrum, are shown in figure 7. 

Another comparison was made for a low Earth orbit (LEO) situation for a representative International 
Space Station (ISS) environment. The electron spectrum and associated calculation were supplied by 
W. Atwell, Technical Research Fellow, Boeing Corporation, Houston, Texas (private communication, 
March 2001), for transport through a semi-infinite aluminum slab (A1 slab) geometry for omnidirectional 
flux subtending 2 jt sr. A Boeing in-house code developed from a series of Monte Carlo calculations was 
used to provide the resultant doses with energy deposition evaluated for an H 2 O medium. Again, the 
comparison is very favorable, as is seen in figure 8. 
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A third comparison calculation was made for geostationary Earth orbit (GEO) conditions based on 
published results of Stassinopoulos (ref. 11) in which the daily omnidirectional integral flux spectrum is 
tabulated along with resultant doses evaluated for energy deposition in an aluminum medium. Figure 9 
shows the results of this comparison. While the electron dose functions compare very favorably, the 
bremsstrahlung curves exhibit considerable differences. Apparently, the larger photon flux at low thick- 
nesses is a consequence of the differential flux evaluation at low energies. Since the energetic photons 
are much more penetrating than electrons, an overestimate of the initial low-energy differential electron 
flux will result in larger bremsstrahlung contributions at greater thicknesses. An investigation of the sen- 
sitivity of results on different evaluations of the derivative was not pursued in this study. 

Sample Calculations and Results 

A description of the program setup, along with examples of optional output, is given here only as an 
illustration of the existing run process and does not pertain to any specific application. The program re- 
quires an input file specifying an initial energy grid and corresponding electron fluxes, differential in 
energy. A maximum energy (normally the highest energy prescribed for the input spectrum) is also speci- 
fied. For the present example, the fluence spectrum of the STS-63 mission used previously has been cho- 
sen, consisting of 29 spectral points with maximum energy of 6 MeV. (Note that the energy range of 
validity for the code in general is 0.01 to 100 MeV.) In the course of the electron transport computational 
procedure, the energy grid points for the electrons are redistributed in accordance with energy loss, with 
density of grid points being greatest at lower energies to better define the slowing down and stopping 
processes. The original number of grid points remains constant. The photon energy grid is automatically 
established for the same number of grid points and increase by logarithmic increments from 0.01 MeV to 
the initial maximum electron energy, and remains invariant throughout the calculation. As an option, a 
photon spectrum may also be included as an initial condition; in the present example, initial photon flux is 
zero. 

A second input file defining the material layers (compositions and thicknesses) is also required. The 
total thickness of the layer and the number of spatial grid points to be included within the layer must be 
specified. The spatial grid points (normally 15 to 20 values) are automatically distributed logarithmically 
for more accurate treatment of the photon transport processes in the material. Values representing thick- 
nesses are actually expressed in terms of areal density, g/cnr (i.e., density x linear thickness). The mate- 
rial composition must be specified in terms of a distinct number of atomic species (presently 1 to 6 ) with 
their respective atomic numbers, atomic weights, and number of atoms of each type in the composite 
molecular formula. These data are then used in the algorithm to compute molecular weights and atomic 
densities per unit mass of material. Data for each layer appear sequentially in the material data file; there 
is presently no limitation on the number of different material layers to be included in a particular calcula- 
tion. A third and final requisite input file contains all of the atomic data for electron stopping powers, 
photon production, and attenuation cross sections, and is only changed for updates to basic atomic prop- 
erties data. 

For the present sample calculation the material slab is composed of three layers: a lightweight (Lt), 
lower molecular weight layer; a heavier (Hv), high molecular weight layer; and a final water layer. A 
polyurethane, (C 5 H 5 NO 2 V, has been chosen as the light material, and a hypothetical tungsten-fiberglass 
composite of assumed composition (WSiCb),, represents the heavy layer. (The assumed formula for the 
composite is based on approximate mass fractions of 0.75 for W and 0.25 for SiOo.) The water layer is 
taken as representative of human tissue. Each layer is specified as having areal density of 0.5 g/cm 2 . 
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The electron stopping powers, associated ranges, and photon production and attenuation coefficients 
are calculated for each material by weighted averaging of the atomic coefficients with respect to the cor- 
responding atomic densities. Figure 10 shows the computed stopping powers and ranges for the three 
selected materials. The photon total attenuation and energy absorption coefficients are given in figure 11, 
and the photon production cross sections are shown in figure 12. 

Two transport calculations were made for the three-layer combination in which the sequence of light 
and heavy materials in front of the tissue layer were interchanged. The calculated electron spectra at the 
material interfaces are shown for the two cases in figure 13. The difference in flux values at the first 
interface is most notable, whereas the flux incident on the final water layer is very similar. The figure 
also illustrates the varying electron energy grid as indicated by the location of symbols along the curves. 
Photon source terms are plotted in figure 14, where photon production at each interface is shown. As 
expected, photon production decreases as electron energies are degraded or as electrons are stopped. 
Most notable is the substantially greater photon production in the heavy layer. Corresponding photon 
fluence spectra at the two interior interfaces are given in figure 15. Largest differences are observed at 
the lower energies (10 to 100 keV); when the layer sequence is such that the middle ply is the lighter 
polyurethane, attenuation at higher photon energies is scarcely noticeable for the thickness specified. 

The variation of dose (energy deposition) in the three-layer slab combinations is presented in fig- 
ure 16. It was specified that 20 grid points be used for the first two strata, while only 5 were used for the 
simulated tissue layer. Again, plot symbols indicate the automatic decadal spacing of the points. The 
very prominent bremsstrahlung dose in the heavy layer is most evident. When the heavy layer is sub- 
jected to the unattenuated environment electrons, relatively more electrons are produced and result in a 
higher dose in the final tissue layer. For this calculation, the corresponding electron dose variations in the 
final layer are quite similar for the two cases. However, other calculations for material layers having 
greater differences in molecular weights have shown larger disparities in final electron tissue dose. 

Concluding Remarks 

The major advantages of the electron procedure here described are its computational speed along with 
its versatility with respect to description of transport in arbitrary materials and/or combinations of materi- 
als. The capability of generating optional electron and photon energy spectra at preselected spatial loca- 
tions is another positive feature that would facilitate numerical diagnostic analysis, and allow for special 
applications where spectral definition is of importance. In its present form, the nominal energy range of 
validity is 0.01 to 100 MeV, but the user should recognize that the experimental data on which the origi- 
nal formulation was based dealt predominantly with measurements for energies <30 MeV and for materi- 
als composed of constituents of intermediate to heavy atomic masses. The rationale for extension to 
lower atomic weights and somewhat higher energies has been discussed; experimental verification in 
these domains would be especially helpful, particularly for ultralightweight polymers designed for space 
flight applications. 

The current version of the code is considered well-suited for calculations of electron energy deposition 
for LEO trapped electrons. For applications involving high-energy beams or the more intense Jovian 
trapped environment, it would be desirable to incorporate pair production and certain photonuclear proc- 
esses in the procedure. Such enhancements would necessarily require extension of the existing database, 
and appropriate modification of the computational grid structure. These improvements are presently 
being formulated for inclusion in future versions of the procedure. 
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Figure 1. CSDA stopping powers for selected elemental species. 
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Figure 3. Ratios of parameterized ranges (Rparam) to CSDA range for selected elemental species. Solid curves are 
present study curve fittings. 
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(a) Carbon (Z = 6). 
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(b) Aluminum (Z = 13). 

Figure 4. Scaled photon production differential cross sections versus electron energy for parameter values of energy 
ratio of emitted photon to incident electron (E v /W). 
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(c) Copper (Z = 29). 



(d) Sliver (Z = 47). 
Figure 4. Continued. 
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(e) Tungsten (Z = 74). 



Electron kinetic energy, MeV' 


(f) Uranium (Z = 92). 
Figure 4. Concluded. 
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Figure 5. Bremsstrahlung total attenuation coefficient for selected elemental species. 
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Figure 6. Bremsstrahlung energy absorption coefficient for selected elemental species. 
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(a) Electron fluence for 10-day mission of STS-63 (394 km at 51.6°). 
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(b) Dose calculation for spectrum in part (a). 

Figure 7. Normal incidence propagation results for electrons in H 2 O. 
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(b) Dose calculation for spectrum in part (a). 

Figure 8. Omnidirectional propagation (2jt sr) for electrons through A1 slab. 
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Electronic kinetic energy, MeV 
(a) Electron flux spectrum for GEO (35 790 km; 70°W). 
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(b) Dose calculation for above spectrum. 

Figure 9. Omnidirectional propagation (2n sr) for GEO electrons through A1 slab. 
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(a) Stopping power. 
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(b) Range. 

Figure 10. Computed stopping powers and ranges for three materials used in sample calculations. 
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Figure 1 1 . Computed photon absorption coefficients for three materials used in sample calculations. 
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(a) Polyurethane. 



(b) WSiC >2 composite. 

Figure 12. Computed photon production cross sections for materials used in sample calculations. 
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Figure 12. Concluded. 
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(b) Sequence WSi02 composite (Hv)/polyurethane (Lt)/water. 

Figure 13. Sample calculation of electron fluence spectra for three-layer slab combinations. 
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(a) Sequence polyurethane (Lt)/WSiC >2 composite (Hv)/water. 



Hv (t = 0.0) 
Lt (t = 0.5) 
H 2 0 (t = 1.0) 


Photon energy, MeV 


(b) Sequence WSi02 composite (Hv)/polyurethane (Lt)/water. 


Figure 14. Sample calculation of photon source spectra for three-layer slab combinations. 
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— ♦ — Lt/Hv 
HV/H2O 



Hv/Lt 

Lt/H 2 0 


(b) Sequence WSiC >2 composite (Hv)/polyurethane (Lt)/water. 

Figure 15. Sample calculation of photon energy spectra for three-layer slab combinations at interior interfaces. 
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(b) Sequence WSiC >2 composite/polyurethane/water. 


Figure 16. Variation of dose (energy deposition) with depth for three-layer slab combinations. 
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